Genetic changes from type I interferons and JAK inhibitors: clues to drivers of juvenile dermatomyositis.

OBJECTIVE: To better understand the pathogenesis of juvenile dermatomyositis (JDM), we examined the effect of the cytokines type I interferons (IFN I) and JAK inhibitor drugs (JAKi) on gene expression in bioengineered pediatric skeletal muscle. METHODS: Myoblasts from 3 healthy pediatric donors were used to create three-dimensional skeletal muscle units termed myobundles. Myobundles were treated with IFN I, either IFNα or IFNß. A subset of IFNß-exposed myobundles was treated with JAKi tofacitinib or baricitinib. RNA sequencing analysis was performed on all myobundles. RESULTS: Seventy-six myobundles were analyzed. Principal component analysis showed donor-specific clusters of gene expression across IFNα and IFNß-exposed myobundles in a dose-dependent manner. Both cytokines upregulated interferon response and proinflammatory genes; however, IFNß led to more significant upregulation. Key downregulated pathways involved oxidative phosphorylation, fatty acid metabolism and myogenesis genes. Addition of tofacitinib or baricitinib moderated the gene expression induced by IFNß, with partial reversal of upregulated inflammatory and downregulated myogenesis pathways. Baricitinib altered genetic profiles more than tofacitinib. CONCLUSION: IFNß leads to more pro-inflammatory gene upregulation than IFNα, correlating to greater decrease in contractile protein gene expression and reduced contractile force. JAK inhibitors, baricitinib more so than tofacitinib, partially reverse IFN I-induced genetic changes. Increased IFN I exposure in healthy bioengineered skeletal muscle leads to IFN-inducible gene expression, inflammatory pathway enrichment, and myogenesis gene downregulation, consistent with what is observed in JDM.

Evaluating Chromatin Accessibility Differences Across Multiple Primate Species Using a Joint Modeling Approach.

Changes in transcriptional regulation are thought to be a major contributor to the evolution of phenotypic traits, but the contribution of changes in chromatin accessibility to the evolution of gene expression remains almost entirely unknown. To address this important gap in knowledge, we developed a new method to identify DNase I Hypersensitive (DHS) sites with differential chromatin accessibility between species using a joint modeling approach. Our method overcomes several limitations inherent to conventional threshold-based pairwise comparisons that become increasingly apparent as the number of species analyzed rises. Our approach employs a single quantitative test which is more sensitive than existing pairwise methods. To illustrate, we applied our joint approach to DHS sites in fibroblast cells from five primates (human, chimpanzee, gorilla, orangutan, and rhesus macaque). We identified 89,744 DHS sites, of which 41% are identified as differential between species using the joint model compared with 33% using the conventional pairwise approach. The joint model provides a principled approach to distinguishing single from multiple chromatin accessibility changes among species. We found that nondifferential DHS sites are enriched for nucleotide conservation. Differential DHS sites with decreased chromatin accessibility relative to rhesus macaque occur more commonly near transcription start sites (TSS), while those with increased chromatin accessibility occur more commonly distal to TSS. Further, differential DHS sites near TSS are less cell type-specific than more distal regulatory elements. Taken together, these results point to distinct classes of DHS sites, each with distinct characteristics of selection, genomic location, and cell type specificity.

Comparative Analyses of Chromatin Landscape in White Adipose Tissue Suggest Humans May Have Less Beigeing Potential than Other Primates.

Humans carry a much larger percentage of body fat than other primates. Despite the central role of adipose tissue in metabolism, little is known about the evolution of white adipose tissue in primates. Phenotypic divergence is often caused by genetic divergence in cis-regulatory regions. We examined the cis-regulatory landscape of fat during human origins by performing comparative analyses of chromatin accessibility in human and chimpanzee adipose tissue using rhesus macaque as an outgroup. We find that many regions that have decreased accessibility in humans are enriched for promoter and enhancer sequences, are depleted for signatures of negative selection, are located near genes involved with lipid metabolism, and contain a short sequence motif involved in the beigeing of fat, the process in which lipid-storing white adipocytes are transdifferentiated into thermogenic beige adipocytes. The collective closing of many putative regulatory regions associated with beigeing of fat suggests a mechanism that increases body fat in humans.

A High-Throughput Mutational Scan of an Intrinsically Disordered Acidic Transcriptional Activation Domain.

Transcriptional activation domains are essential for gene regulation, but their intrinsic disorder and low primary sequence conservation have made it difficult to identify the amino acid composition features that underlie their activity. Here, we describe a rational mutagenesis scheme that deconvolves the function of four activation domain sequence features-acidity, hydrophobicity, intrinsic disorder, and short linear motifs-by quantifying the activity of thousands of variants in vivo and simulating their conformational ensembles using an all-atom Monte Carlo approach. Our results with a canonical activation domain from the Saccharomyces cerevisiae transcription factor Gcn4 reconcile existing observations into a unified model of its function: the intrinsic disorder and acidic residues keep two hydrophobic motifs from driving collapse. Instead, the most-active variants keep their aromatic residues exposed to the solvent. Our results illustrate how the function of intrinsically disordered proteins can be revealed by high-throughput rational mutagenesis.

Causal Genetic Variation Underlying Metabolome Differences.

An ongoing challenge in biology is to predict the phenotypes of individuals from their genotypes. Genetic variants that cause disease often change an individual's total metabolite profile, or metabolome. In light of our extensive knowledge of metabolic pathways, genetic variants that alter the metabolome may help predict novel phenotypes. To link genetic variants to changes in the metabolome, we studied natural variation in the yeast Saccharomyces cerevisiae We used an untargeted mass spectrometry method to identify dozens of metabolite Quantitative Trait Loci (mQTL), genomic regions containing genetic variation that control differences in metabolite levels between individuals. We mapped differences in urea cycle metabolites to genetic variation in specific genes known to regulate amino acid biosynthesis. Our functional assays reveal that genetic variation in two genes, AUA1 and ARG81, cause the differences in the abundance of several urea cycle metabolites. Based on knowledge of the urea cycle, we predicted and then validated a new phenotype: sensitivity to a particular class of amino acid isomers. Our results are a proof-of-concept that untargeted mass spectrometry can reveal links between natural genetic variants and metabolome diversity. The interpretability of our results demonstrates the promise of using genetic variants underlying natural differences in the metabolome to predict novel phenotypes from genotype.

Heterochronic meiotic misexpression in an interspecific yeast hybrid.

Regulatory changes rapidly accumulate between species, and interspecific hybrids often misexpress genes. Hybrid misexpression, expression levels outside the range of both parental species, can result from cis- and trans-acting regulatory changes that interact abnormally in hybrids. Thus, misexpressed genes may contribute to hybrid sterility. However, in the context of a whole organism, misexpression may not result directly from cis-trans interactions but rather indirectly from differences between hybrid and parental abundance of cell types. Here we eliminate the confounding effects of cell types by examining gene expression in a sterile interspecific yeast hybrid during meiosis. We investigated gene expression of the yeasts Saccharomyces cerevisiae, S. paradoxus, and their hybrid at multiple meiotic stages. Although the hybrid and parents exhibit similar changes in expression levels across meiosis, the hybrid meiotic program occurs earlier than either parent. The timing change produces a heterochronic pattern of misexpression during midmeiosis. Coincident with the timing of misexpression, we find a transition from predominantly trans-acting to cis-acting expression divergence and an increase in the number of opposing cis-trans changes. However, we find no direct relationship between opposing cis-trans changes and misexpression. Contrary to the notion that cis-trans interactions cause misexpression, a heterochronic shift in the normal meiotic gene expression program produces patterns of misexpression in an yeast hybrid. Our results imply that temporal dynamics of single cell types is important to understanding hybrid misexpression and its relationship to cis-trans interactions.